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Scientific  Progress 


The  scientific  progress  will  be  reported  by  subtask  as  contained  in  our  original 
proposal.  In  mid-August  we  gave  a  detailed  briefing  to  Dr.  Kenneth  Dodd  of  WRAIR  on 
the  technical  accompishments  of  the  project  to  date  and  to  Col.  Sedge,  RAD-3,  on  the 
technical  and  programmatic  progress.  The  materials  presented  in  those  briefings  are 
attached. 
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Research  Topic  1:  Pathology  of  Lung  Injury 

Task  1.  Construction  of  a  database  for  lung  injury  pathology 

Subtask  1.1  Database  specification.  Working  together  with  the  WRAIR  COR,  we 
have  formulated  a  specification  of  the  entries  to  be  included  in  the  pathology  database. 
That  specification  was  submitted  with  the  1st  quarterly  progress  report.  The  computer 
program  for  entering  new  pathology  data  will  follow  this  format. 

The  existing  Albuquerque  pathology  data  collected  and  analyzed  under  this  Task, 
however,  does  not  include  this  level  of  detail.  Consequently,  this  data  must  be 
characterized  differently  to  avoid  introducing  inordinate  speculation.  A  compromise  was 
arrived  at  in  which  a  new  entry  is  included  for  each  anatomical  location  (larynx,  trachea, 
five  lung  lobes,  and  13  sections  of  the  gastrointestinal  tract)  that  indicates  injury  on  a 
yes/no  basis.  This  will  allow  the  existing  data  to  be  correlated  and  will  not  conf  ict  with 
more  complete  data  entered  later. 

After  all  of  the  pictorial  data  on  lung  injury  has  been  analyzed,  we  will  work  together 
with  the  WRAIR  COR  to  refine  the  injury  description  in  each  of  the  lung  lobes. 

This  subtask  is  now  complete. 

Subtask  12  Software  development.  The  pathology  database  software  consists  of  two 
parts.  The  first  is  a  relational  database  for  the  entries  in  the  specification  sheet  and  related 
data  on  the  blast  exposure  (text  data).  The  second  is  a  program  for  entering,  retrieving, 
and  manipulating  graphical  representations  of  the  lung  necropsy  pictures. 

The  software  to  enter,  modify,  and  browse  the  text  data  associated  with  the 
Albuquerque  tests  has  been  written  in  FORTRAN  and  C  Languages  using  the 
commercially  available  package  of  database  routines  made  by  INFORMIX.  The  database 
now  includes  all  of  the  blast  and  pathology  entries  necessary  to  store  the  existing  data.  It 
was  necessary  to  add  database  records  to  store  lengthy  comments  that  were  part  of  the 
Albuquerque  posting  sheets. 

This  part  of  the  subtask  is  now  complete. 

Future  animal  tests  will  report  necropsy  data  according  to  the  newly  designed,  more 
quantitative  pathology  sheets.  We  are  modilying  the  database  software  to  allow  the  user  to 
enter,  modify,  and  browse  through  data  in  this  form.  We  understand  that  WRAIR  would  like  to 
re-score  all  of  the  old  data  into  the  new  format,  so  that  there  might  eventually  be  a  single, 
consistent  data  format. 
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Subtask  1.3  Compilation  of  data.  All  of  the  data,  magnetic  tape,  necropsy  scoring 
sheets,  strip  chart  records,  field  notebook  records,  photographs,  and  written  reports  that 
are  available  to  us  have  been  compiled  and  cataloged.  There  are  two  areas  where  data  is 
incomplete:  digitized  pressure  traces  from  the  "double  peak"  study  and  the  original 
pathology  sheets  from  one  of  the  Albuquerque  tests.  The  blast  data  requires  the  tape 
recorder  used  at  the  Kirtland  Test  Site,  while  EG&G  has  the  pathology  data.  WRAIR  has 
requested  both  from  the  operators  of  the  Test  Site.  Since  the  data  has  not  been  made 
available  soon,  we  will  not  be  including  it  in  this  version  of  the  database. 

The  pathology  data  has  been  organized  by  test  series  and  by  animal  within  test  series. 
All  photographs  have  been  sequentially  numbered  so  they  can  be  cross-referc  ':ed  within 
the  database.  We  have  made  high-quality  color  reproductions  of  all  of  the  data  and 
returned  the  original  data  to  WRAIR  for  safekeeping. 

The  data  comprises  tests  of  over  600  animals  and  the  hardcopy  form  consists  of  over 
5000  pages  of  text  and  pictures. 

This  subtask  is  complete. 

Subtask  1.4  Qualification  of  the  data.  All  of  the  data  has  been  reviewed  for 
consistency  and  completeness.  No  obvious  inconsistencies  have  been  found,  although 
anomalies  have  been  seen,  for  example,  one  animal  severely  injured  when  five  others  with 
the  same  exposure  are  not  injured.  These  differences  may  reflect  animal  variability  or  may 
later  be  found  to  be  a  reporting  error.  Resolution  of  these  questions  must  await  the 
findings  of  the  model  tasks. 

The  data  is  satisfactorily  complete.  The  earlier  tests  were  less  detailed  in  reporting 
pathology,  but  that  was  expected,  since  it  was  not  a  requirement  at  the  time.  The  pathology 
data  that  EG&G  has  refused  to  release  is  an  obvious  problem,  but  it  affects  the  graphical 
aspects  rather  than  the  text  entries. 

This  subtask  is  complete. 

Subtask  1.3  Data  entry.  All  of  the  pathology  text  data  that  has  been  made 
available  to  us  has  been  entered. 

This  part  of  the  subtask  is  complete. 
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We  will  work  with  the  WRAIR  COR  during  the  next  quarter  to  determine  the  needs  for 
electronic  form  of  the  graphical  and  photographic  data. 

Subtask  1.6  Reporting  and  distribution.  Delivery  of  these  items  is  not  planned  for 
another  year,  however,  the  database  has  already  answered  some  immediate  needs.  At  the 
request  of  the  WRAIR  COR  we  analyzed  the  existing  entries  to  determine  the  probable 
validity  of  the  Bowen  correlation  in  predicting  injury.  The  preliminary  results  show  that 
injury  occurs  at  blast  levels  a  factor  of  two  or  three  lower  than  that  predicted  by  the  Bowen 
curve  for  threshold  lung  injury. 

Subtask  1.7  Database  maintenance.  Work  has  not  begun  on  this  subtask. 
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Task  2.  Multi-dimensional  model  of  thorax  response 

Subtask  2.1  Multi-dimensional  free-fleid  load  description.  Load  data  taken  in  the 
free-field  from  Lambdroid  and  chest  wall  mounted  gages  and  well  as  previous  EITACC 
calculational  results  have  been  compiled.  The  results  support  the  following  correlation.  The 
blast-side  loading  is  characterized  by  the  reflected  wave,  while  the  leeward  side  is  comparable  to 
the  incident  wave  and  is  delayed  in  time  by  about  1  to  2  msec.  Further  refinement  of  this 
description  is  not  justified  by  the  current  data.  Since  there  are  no  new  experiments  being 
conducted  at  Albuquerque,  it  would  be  deshable  to  make  more  EITACC  CFD  calculations  in 
the  future. 

This  subtask  is  complete. 

Subtask  2.2  Multi-dimensional  thorax/lung  model.  Previous  Finite  Element 
Modeling  (FEM)  of  the  thorax/lung  has  modeled  both  the  chest  wall  and  the  lung  with 
elements,  which  led  to  aphysical  solutions  because  of  the  large  density  difference  between 
the  chest  wall  and  the  lung  parenchyma.  A  one-dimensional  model  was  developed  that 
represented  the  chest  wall  mass  as  a  single  node,  resolving  the  previous  problems. 

The  technique  of  treating  the  chest  wall,  and  other  solid  boundaries  of  the  lung,  as 
part  of  the  nodes  has  been  implemented  in  a  two-dimensional  FEM.  The  resulting 
solutions  are  well  behaved  and  agree  with  the  localized  one-dimensional  results. 

We  have  ordered  and  received  an  anatomical  model  of  the  sheep  lung  to  provide  specific 
guidance  for  refinement  of  the  FEM.  Based  on  that  model,  a  new  FEM  grid  has  been  selected 
that  better  resolves  the  stresses  at  the  pleural  surface  and  whose  elements  have  an  initial  shape 
with  less  distortion.  Both  of  these  changes  enhance  the  calculational  accuracy.  Further 
refinements  of  the  chest  wall  treatment  are  being  implemented  to  allow  proper  description  of 
the  rigid  nature  of  the  thoracic  cage,  without  compromising  the  tangential  motion  of  the  pleural 
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Task  3.  Validation  against  the  pathology  database 

Subtask  3.1  Generalization  of  parenchymal  work  correlate.  Work  on  this  subtask 
has  not  begun. 

Subtask  3.2  Validation  of  overall  injury  characterization.  Work  has  begun  to 
validate  the  INJURY  code  against  the  overall  injury  scores  in  the  pathology  databtise.  The  fi''st 
step  is  to  modify  the  program  to  read  the  pressure  data  in  .JIF  format,  the  new  standard  for 
exchanging  data  between  the  software  modules.  In  the  process  of  making  this  change,  it  was 
clear  we  needed  to  rearrange  the  functional  aspects  of  the  program  to  allow  ready  access  to  the 
model  parameters  and  to  provide  better  output  formats  for  analyzing  the  data. 

The  numerical  solution  algorithm  was  re-examined.  It  was  discovered  that  the  5th  order 
Runge-Kutta  method,  which  is  ideal  for  most  differential  equations,  is  ill-suited  to  the  sharply 
peaked  pressure  traces  associated  with  blast  waves.  We  experimented  with  an  exact  analytical 
solution  method,  but  in  the  end  found  that  a  fully-implicit,  time-averaged  scheme  was  not  only 
just  as  accurate,  but  was  more  easily  implemented.  The  code  now  computes  the  injury  models 
in  a  fraction  of  the  previous  time  (for  the  same  accuracy). 

We  have  begun  to  link  the  pathology  and  blast  databases.  The  INJURY  program  requires 
a  pressure-time  history  to  make  a  calculation  of  the  body  dynamics.  By  linking  the  databases 
we  are  linking  the  observed  injury  with  the  measured  pressure  record.  A  single  unified  database 
is  beyond  the  scope  of  the  current  effort,  because  major  reorganization  of  the  VU  file  format 
would  be  required.  Instead,  we  will  create  separate  .JIF  files  for  each  shot  that  can  be  used 
with  INJURY. 

Subtask  3.3  Validation  against  pathology  distribution  data.  Work  on  this  subtask 
has  not  begun. 

Subtask  3.4  Model  refinement.  Work  on  this  subtask  has  not  yet  begun. 
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Task  4.  Development  of  a  human  thorax/lung  model 

Subtask  4.1  FEM  human  thorax/lung  model.  Work  this  subtask  has  not  yet 
begun. 

Subtask  4.2  Modifled  injuiy  predictions.  Work  on  this  subtask  has  not  yet  begun. 

Subtask  4.3  Human  model  into  INJURY  software.  Work  on  this  subtask  has  not  yet 
begun. 

Task  5.  Extension  of  Pathology  Database  to  GI  Injury 


Subtask  5.1  Data  Entry.  The  database  was  extended  to  include  fields  for  the 
gastrointestinal  tract.  All  of  the  available  Albuquerque  text  data  has  been  entered. 

This  subtask  is  complete. 

Subtask  5.2  Correlation  of  Injury.  Work  has  begun  to  correlate  the  gastrointestinal 
injury  with  blast  wave  parameters.  Preliminary  results  will  be  passed  to  the  Defense  Nuclear 
Agency  as  soon  as  possible. 

Subtask  5.3  Reporting  and  Distribution.  Work  is  not  planned  untH  next  year. 
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Research  Topic  2:  Critical  Parameters  in  Complex  Waves 

Task  1.  Complete  Description  of  Complex  Waves 


Subtask  1.1  Compilation  of  load  data  in  enclosures.  A  review  of  the  available  load 
data  for  blasts  within  enclosures  has  been  performed.  The  only  data  sufficiently  detailed 
to  serve  the  purposes  of  injury  modeling  is  that  collected  at  the  Albuquerque  Test  Site  for 
WRAIR.  The  Navy  charge/volume  scaling  rules  are  more  applicable  to  the 
static  overpressure  in  unvented  spaces.  The  information  pertinent  to  the  use  of  this  data 
for  load  determination  has  been  entered  into  the  blast  database  and,  where  possible,  the 
traces  have  been  digitized  for  future  analysis. 

This  subtask  has  been  completed. 

Subtask  1.2  CFD  calculations.  JAYCOR’s  EITACC  computer  program  has  been 
applied  to  the  simulation  of  single  charge  explosions  within  an  enclosure.  The  results  have 
been  compared  with  the  measurements  from  the  Albuquerque  Test  Site  to  determine  the 
range  of  validity  of  the  model  and  with  the  method  of  images  computer  program 
COMPLX.  The  results  show  that  CFD  methods  can  be  used  to  predict  much  of  the  detail 
of  the  complex  waves  and  to  provide  an  analytical  source  of  data  for  loading  that  can  be 
used  in  injury  prediction  models.  The  details  of  this  investigation  have  been  described  in 
a  draft  technical  report  delivered  in  January,  1990. 

This  work  will  eventually  be  incorporated  into  a  single  report  on  complex  waves  that 
includes  the  results  of  Task  2. 

This  subtask  has  been  completed. 

Subtask  1.3  Critical  parameters  of  blast.  Work  has  not  begun  on  this  subtask. 

The  CFD  calculations  can  be  quite  time  consuming.  During  the  coming  quarter, 
JAYCOR  is  taking  delivery  on  new,  high-speed  workstations  that  are  almost  a  factor  of  ten 
faster  than  our  current  workstations.  We  plan  to  defer  this  subtask  (and  sub  task  1.4)  until  the 
following  quarter,  when  the  hardware  will  be  fully  operational  and  the  calculations  can  be 
made  more  cost  effectively.  The  delay  will  also  allow  us  more  time  to  analyze  the  existing  data 
and  the  predictions  of  the  COMPLX  code. 

Subtask  1.4  Load  distribution  for  complex  waves. 

Work  has  not  begun  on  this  subtask.  See  previous  subtask. 
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Task  2.  Incorporate  findings  into  COMPLX 

Subtask  2.1  Extend  the  COMPLX  model.  Work  has  begun  to  analyze  the  predictions 
of  the  COMPLX  code.  Several  areas  of  improvement  have  been  identified.  First,  the 
current  algorithm  for  determining  the  intersection  of  a  blast  ray  with  a  reflecting  surface  was 
found  to  be  sensitive  to  numerical  precision.  A  new  algorithm  has  been  derived  and  is  being 
tested.  Second,  the  previous  algorithm  terminated  the  image-generating  cycle  after  a  fix 
number  of  reflections.  Further  analysis  revealed  that  many  higher-order  reflections  are 
important  to  the  pressure  signal  in  the  10  to  20  msec,  time  frame.  It  appears  that  part  of  the 
discrepancy  with  data,  that  had  earlier  been  attributed  to  pressurization  effects,  are  actually 
explained  by  multiple  reflections. 

Subtask  2.2  Validation  against  CFD  model.  Work  has  begun  to  compare  the 
COMPLX  calculations  with  the  EITACC  CFD  calculations.  The  agreement  is  qualitatively 
good,  but  there  are  clearly  features  that  only  the  CFD  calculations  can  capture.  The 
investigation  will  focus  on  the  physical  origin  of  those  differences  and  whether  they  can  be 
incorporated  into  the  BWAVES  code. 


Contract  DAMD17-89-C-9150 


Midterm  Peport 


Task  3.  Critical  parameters  of  injury 
Subtask  3.1  Evaluate  multi-dimensional  aspect  of  injuty. 
Subtask  32  Validation  of  the  extended  COMPLX  iiyury  model. 
Subtask  3.3  Determination  of  critical  parameters. 

Subtask  3.4  Upgrade  INJURY  software. 


Work  on  this  Task  has  not  begun. 
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August  16, 1990 


BASIS  OF  CURRENT  ARMY  DOCTRINE 


Non-auclitory  limit  (Mil  Std  1474B  at  5  msec)  Ps  =  0.34  atm 
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No  guidance  on  application  to  complex  waves 
[Pb  does  not  correlate] 
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Findings:  general,  but  unvalidated 
Products:  reports,  user-oriented  software 


GOALS  OF  CURRENT  CONTRACT 
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Extend  pathology  database  to  Gl 

Evaluate  importance  of  jDjast  on  iDattlefield 

Estimate  the  medical  and  performance  consequences 


Task  1.  Lung  Pathology 

ORGANIZE  ALBQ  DATA  INTO  USABLE  FORM 
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Data  Qualification 

Blast  parameters 
Pathology  re-grading 

Product:  Extension  of  Baker’s  relations  for  Mach  region 
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Figure  2.  Comparison  of  measured  peak  incident  pressure 
with  value  calculated  from  a  modification  of  Baker's  relations 
that  accounts  for  wave  interaction  in  the  Mach  stem. 
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Figure  3.  Comparison  of  mi 
with  values  calculated  from 


Task  1.  Lung  Pathology 
QUANTIFY  TRENDS 
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Products:  Quick-Look  Report 

Statistical  analysis  in  INJURY 


Animal  Weight  Variation 


Figure  1.  Variation  of  animal  weight  with  animal  number 
Gaps  represent  animals  not  used  in  subset.  Since  animal 
numbers  are  assigned  sequentially,  a  decrease  in  weight  in 
later  tests  can  be  seen. 
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Figure  5a.  Comparison  of  observed  occurrence  of  no  lung 
injury  with  analytical  function. 


Injury  Frequency  Correlation 

Slight  Lung  Injury 
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Figure  5b.  Comparison  of  observed  occurrence  of  slight 
lung  injury  with  analytical  fuction. 
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Figure  5c.  Comparison  of  observed  occurrence  of  moderate 
lung  injury  with  anlaytical  function. 
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Figure  5d.  Comparison  of  observed  occurrence  of  severe 
lung  injury  with  analytical  function. 
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Figure  5e.  Percent  incidence  of  lung  injury  as  a  function 
of  the  Bowen  Pressure  as  correlated  with  analytical  functions. 


Injury  Occurrence  Correlation 

Tacheal  Injury 
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Figure  4e.  Percent  incidence  of  tracheal  injury  as  a  function 
of  the  Rowen  pressure  as  correlated  with  analytic  functions. 
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Figure  6e.  Percent  incidence  of  gastrointestinal  injury 
as  a  function  of  Bowen  pressure  as  correlated  with 
analytical  functions 


Baker  friedL  File  Quit 

Ca leu Ici'te -and  Show  results 


'•/  0 


in 

C9 

(S 

IS 

. 

(S 

IS 

o 

(4 

CD 

(S 

CD 

in 

in 

. 

C9 

s 

(S 

in 

• 

in 

(0 

•t* 

'  ■ 

• 

■ 

H 

cn 

00 

. 

01 

E 

rt! 

II 

II 

II 

II 

II 

II 

II 

II 

H 

it! 

V) 

a; 

O* 

H 

cm 

M 

ii 

± 

D 

« 

s 

•P 

o 

M 

03 

D 

M 

in 

in 

M 

Pi4 

05 

05 

in 

K 

H 

> 

1*4 

05 

'h5 

Or 

3 

H 

o 

U 

z 

O 

$4 

n 

K 

o 

Oh 

Oh 

M 

M 

H 

Oh 

IH 

Z. 

1 

y 

y 

X 

H 

H 

HH 

c 

1 

« 

o 

O 

X 

Oi 

1 

c: 

M 

Z 

<£ 

oi; 

3 

1 

X 

X 

M 

<c 

M 

3 

a 

0 

1 

0 

o 

X 

Oh 

1=1 

z 

n 

9 


m 


n 


o 

—  • 

■■  CO 

TTT 

19 

jK 

- 

1 

“ 

1 

nj* 

_ 

1 

V 

H 

:  ® 

s 

H 

■ 

D 

U 

Ga 

•H 

D 

in 

H 

C 

(0 

0 

t. 

H 

Ol 

■9 

-o- 


1iriiiiili]iiiiniMi|iiiiiiiii|iiiii.iiii 

-0  ’eST  0  ’Oi  •  *  0 '0T- 


31  W 
H  Vl  ,  H 
fl) 


-fc-i- 


36 


Calculate  an.d  Sllow  results  Special  Messagearea 


Task  1.  Lung  Pathology 
VALIDATE  BIOMECHANICAL  MODELS 
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Work  correlate  extension,  new  INJURY 


Figure  4-1 .  APC  specification  and  sensor  location. 
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Figure  4-2.  Compu.  -tlonal  cell  network  for  calculating  the  blast  in  an  ARC. 
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Figure  4-5.  Comparison  between  calculated  and  measured  pressure  histories  from  a  blast  of 
1  lb  of  C-4  in  an  APC. 
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(c)  In  crew  compartment 

Figure  5- 1 .  Pressure  time  history  from  a  blast  of  1  lb  of  C-4  in  an  APC. 
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-igure  8-1 .  Definition  of  free  field  sensor  locations  in  the  field  bunker. 
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Figure  8-2.  COMPLX  results  for  front  sensor,  1  lb  C4  in  bunker. 
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Product:  Topical  report,  upgraded  BWAVES 


Affects  pressurization  and  reflection  at  later  time 
Changes  effective  propagation  speed  (temperature) 


TASK  3.  COMBINED  INJURY 
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Will  incorporating  combined  injury  effects  increase  this  distance  by  at  least  100m? 
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COMBINED  INJURY  MODELS 
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Mechanistic  Models 


MECHANISTIC  MODEL 
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SELECTION  OF  PARAMETERS 
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SELECTION  OF  PARAMETERS  (Cont’d) 
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1000  rad  Bq  =  0.75  T,  =  200  hrs 

I  ,1  . 


CALIBRATION  AGAINST  RADIATION 
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DELIVERABLES 


(n 

(0 
«  aBis 

0 

03 

£ 

■D 

0 

0  .N 

0  C 
CC  0 

•9  03 

O 
0 


0 

0 

$ 

4:; 

o 

0 

o 

*0 

□ 

C 

0 

E 


0 
0 

Q  cf 

>. 

03 
O 

O 
sz 

4-» 

0 
CL 


0  .h: 


0 

Ll 

I 

0 

0 

(f)  vL 
03  U- 
C 


0 

■Q 

O 

E 

X 

0 

o 

sz 


0 

0 

D 


•O 

c 

Ll 


o 


03 

.E 

*0 

•o 

o 


13 

'c 


}- 

cn 

3 


0 
■o 
o 

E 

X 

0 

w 

o 

sz 

i « 

0  c 

0  c 
sz  o 
0  sz 


LU  UJ 
LL  LL 


C 

o 

’■+—• 

0 

3 

0 

> 

■0 

■o 

o 

E 

*o 

0 

.N 

"0 


0 

4-» 

Z3 


0 

C 

o 


O 

0  o 
0  < 


13 

w  0  ^ 

0^3 
CL 


O 


0 


>> 


0 

C 

0  QT  ■  \_y  •— 

O  §  LU  O  £ 


«  £  o  #  f 

■o 
0 
C 


X 

0 

CL 

E 

o 

o 


S2 

E 

o 

O 


63 


Blast  importance  on  battlefield 

Biological  response  model  O  working  document 

•  Delivered 


DEm^RABLES 


>» 

OD  CO 
O  LU 

O  > 

a.  cn 


00  ^ 
CO 

UJ  < 

o  1 

O  Q 

<  < 


Q 

®  C 
<  10 
?  E 

Q  O 

<  O 


CjJ 

Q  O 

CL  ^ 

$  < 

^  CL. 


o  Working  document 
•  Delivered 


SCHEDULE 


>^ 

Q) 

H  o 

J  o 

p.  ^ 


Q  cd 

W  a. 

.. 

U  -- 

^  o 


0 

CO 

C  0 
O  CO 
Q.  cd 
0  CD 

CC  -2 


cd  xj 

goo 

O  >^2 

^  o  :□ 
o 

cd  X 

®Q.  2 
0-^0 
c 

•—  f— “ 

"Tri  cd 

C  CD  C 

Q<  P 

0  C  I 
c  o  ^ 
0  ’-P  -L 
F  cd 
.F  -D  ^ 

■O  ^ 
p:  0  Jil 

S> 


CM 

CO 

js::  0 
CO  h- 
0 
H 


66 


SCHEDyLE 


O)  s 
o  ^ 


o  ^ 

Cl  ^ 
O  ^ 

h- 


67 


SCHEDULE 

Research  Topic  2: 
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